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ABSTRACT 

We assess the fraction of the Milky Way’s circumgalactic medium (CGM) eluding detection due to 
its velocity being similar to gas in the disk. This is achieved using synthetic observations of the CGM 
in a simulated MW-mass galaxy that shows similar CGM kinematics to the MW and external L^L* 
galaxies. As viewed by a mock observer at a location similar to the Sun, only 50% (by mass) of the gas 
moves at high velocity (|plsr| >100 km s -1 or |pdev| >50 km s -1 ) in the simulated CGM and would 
be observable. The low velocity gas is thermodynamically similar to the high velocity gas, indicating 
the 50% observable fraction is applicable to spectral lines from the radio to the ultraviolet. We apply 
the observable mass fraction (50%) to current estimates of the MW’s CGM, and find a corrected total 
mass of 2.8xlO 8 M© for gas below 10 6 K within ^ 15 kpc (this excludes the Magellanic System). This 
is less than the total mass of the CGM extending out to ^150 kpc in other L^L* galaxies. However, 
we find similar O VI column densities when the discrepancy in path length between the MW and 
external galaxies is considered. The coherent spatial and kinematic distribution of low velocity gas 
in the simulated CGM suggests that current H I observations of the MW’s CGM may miss large low 
velocity H I complexes. In addition, current mass estimates of the MW’s CGM based on high-velocity 
line observations with distance constraints may miss a non-negligible fraction of gas in the outer halo 
which can be obscured if it moves at a velocity similar to the gas in the lower halo. 

Subject headings: Galaxy: kinematics - Galaxy: halo - methods: numerical 


1. INTRODUCTION 

Studies of the circumgalactic medium (CGM) help us 
to understnad the origin of a galaxy’s star formation 
fuel and the impact of galactic feedback. Assessing the 
amount and state of the baryons in galaxy halos is a 
key active area in galaxy formation research. Due to its 
proximity, the Milky Way (MW)’s CGM has been stud¬ 
ied in the most detail, and this has helped us interpret 
the nature of the multiphase medium and the potential 
interactions with the galactic disk. 

To calculate the mass of the multiphase CGM in the 
MW, numerous spectroscopic observations towards high- 
velocity line features have been conducted. High velocity 
clouds (HVCs), mapped in H I surveys via the 21 cm hy- 
perfine transition line, trace the coldest component of the 
Galactic halo. HVCs have been found to largely reside 
within 15 kpc of the disk with the exception of the com¬ 
plexes associated with the M agellanic S ystem which are 
at least at ^55 kp( 


man et al. 

2012 

non-Magel 
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2012, hereafter PPJ12). The H I mass of the 


and thus is a very small fraction of the Galaxy’s baryons. 
Observations of the ion absorption lines in the spectra 
of background quasars show an order of magnitude more 
mass is associated with t he warmer ionized gas (T^10 4 ~ 6 
K) in the Galactic halo ( Sembach et al.|2003||Shull et al.| 


2009[ PPJ12). X-ray observations of hot gas (TMO 0 K) 


m emission and absorption (O VII, O VIII) do not have 
strong constraints from spectral line observations, but 
currently indicate a mass of ^10 10 M© if extended out 
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to 250 k pc (| And erson fe Bregman||2010l iMiller & Breg- 
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Tiie MW’s CGM can be compared to that of other 
L^L* galaxies now with the increase in absorption line 
probes made available with the Cosmic Origins Spectro¬ 
graph (COS). Discrepancies in both mass and column 
density have been found. At first glance, the total mass of 
the MW’s CGM with T<10 6 K appears to be 1-2 orders 
of magnitude lower than the CGM mass of other 
spirals (~10 10 M, 
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Stocke et al.|2013l Werk et al.|2014 ). 


As for the column density, Tumlinson et al.|(|201ip found 
a typical value of log Novi=14.5 in the CGM of L~L* 
galaxies with specific star formation rate (sSFR)>10 -11 
yr -1 in the COS-Halos samp le, while the MW’s CGM 
instead has log (Novi )=14.0 (|Sembach et al. 2003 J) wi th 


sSFR~2-6xlO 11 yr 1 (SFR^l -4 M© yr 1 


Robitaille fe 

Whit ney|2QlQ||F)iehl et al.|2006| log M* ^10.8, |lVlclVlillan 
201 ip. These differences between the MW and external 
L^L* galaxies raise several questions. Is the MW’s CGM 
an outlier with much less mass than other L^L* spirals? 
Or are our observations of the MW’s CGM still incom¬ 
plete? 

Observations of the MW’s CGM specifically suffer from 
obscuration by the dense foreground disk gas that can 
be moving at a similar radial velocity as the gas in the 
MW’s CGM. The disk gas causes studies of the MW’s 
CGM to apply velocity cutoffs with commonly adopted 
values of IlysrI < 90 — 100 km s" 1 (Wakker |l991a 


Sembach et al.| 2003 Lehner et al. 201!? and PPJ12). 


It is also common to exclude gas with deviation ve¬ 
locity I^devI < 50 — 75 km s -1 which is the amount 
a cloud’ s velo city devi ates from a simple disk rotation 
model (Wakker 1991al PPJ12). However, estimates of 
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the obscured MW's CGM mass due to these velocity cut¬ 
offs are rarely presented since it is usually impossible to 
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distinguish signals of low velocity halo gas from that of 
disk gas. 

In this work, we conduct synthetic observations of 
the CGM of a simulated MW-mass galaxy. Observa¬ 
tional conditions of the MW’s CGM are reproduced so 
that the influence of foreground disk material on the 
observed CGM properties can be evaluated. Following 
the definitions in Wakker (1991a), two velocity cutoffs 
(I^lsrI > 100 km s -i and I^devI > 50 km s -1 ) are 
adopted to define the high velocity gas that is observable 
in the simulated CGM. We do not perform sophisticated 
modelling of the gas ionization state, but rather approxi¬ 
mate the CGM gas phases by several temperature ranges: 
cold, T<10 4 K; warm, TW10 4-5 K; warm-hot, dWIO 5-6 
K; hot, T>10 6 K. These loosely correspond to gas traced 
by atomic hydrogen (H I), low and intermediate ions 
(e.g., Si II, Si III, Si IV, C II, C III, C IV), high ions 
(e.g., O VI), and X-ray emitting/absorbing ions (e.g., 
O VII, O VIII), respectively. With this setup, we can 
understand the CGM distribution in high and low veloc¬ 
ity regimes at various phases, which will shed light upon 
the question of whether the MW’s CGM is abnormal or 
current observations are incomplete. 

This paper is organized as follows: Section [2] briefly 
describes the setup of the numerical simulation and gas 
kinematics of the simulated CGM; Section [3] outlines the 
methods we use to conduct synthetic observations on the 
simulated CGM with a view from within the simulated 
disk. The properties of obscured low velocity gas in the 
simulated CGM through spectroscopic observations are 
shown in Section [4j and comparisons with current ob¬ 
servations of the MW and other L^L* galaxies are dis¬ 
cussed in Section [5] We conclude in Section [6] with our 
main findings. 


the disk in the form of cold clouds (H I). Another source 
of H I is the stripped gas from satellites as they move 
through the simulated CGM. Inflows and the stripped 
satellite gas contribute comparable amounts of H I to 
the halo. Galactic winds also affect the gaseous CGM. 
Large-scale collimated outflows driven by supernova feed¬ 
back have been reported in the simulated galaxy (also see 
Section 3]); these winds generally carry very hot metal- 
enricheagas and move at high velocities extending from 
200 km s- 1 to 1000 km s” 1 (PPJ12). 

3. SYNTHETIC OBSERVATIONS 

In this paper, the simulated disk is defined as a cylin¬ 
drical region with a size of R<18 k pc a nd |z| <2 kpc 


following the terminology used in Joung et al. (2012) 
The simulated CGM defines a spherical volume witn a 
radius of 250 kpc with the central simulated disk manu¬ 
ally masked. We use the yt analysis toolkit (Turk et al 
2011) to study the properties of the simulated CGM from' 


a view within the simulated disk, and perform analyses 
on the CGM gas with all the temperature available in 
the simulation unless otherwise specified. The synthetic 
observations account for the galactic differential rotation 
of the simulated disk and the off-center location of the 
mock observer (Local Standard of Rest; LSR), which is 
the case for observations of the MW’s CGM. 

We construct a mock Galactic coordinate system anal¬ 
ogous to that of the MW, and specify observing sight- 
lines using galactic longitude l and galactic latitude b. 
We then build a disk rotation model by synthesizing H I 
21cm emission from the simulated disk in order to deter¬ 
mine the velocity ranges allowed by disk rotation, within 
which low velocity gas in the simulated CGM at similar 
velocities will be obscured. 


2. SIMULATION OF A MW-MASS GALAXY 

Our investigation of the gas kinematics in the CGM 
is based on a high-resolution cosmological simulation of 
a MW- mass galax y. The simulation was performed by 
Joung et al. (2012) with Enzo , an Eulerian hydrodynam¬ 
ics adaptive refinement code. A low-resolution run was 
first performed which identified four MW-mass halos in 
Local-Group-like volumes; then one of the halos was re¬ 
simulated, achieving a mass resolution of 10 5 M 0 and 
spatial resolution of 136-272 pc physical (or better at all 
times). The code includes metallicity-dependent cooling, 
photoionization by UV extragalactic radiation, shield¬ 
ing of UV radiation by H I, and photoelectric heating 
of the interstellar medium. Star formation and stellar 
feedback (but no AGN feedback) are implemented. We 
refe r the reader to th e work done by Joung et al. (2012) 
and Fernandez et al.| (|2012|) for detailed information on 
the performance of the simulation, and only report their 
key results here. 

The total mass of the simulated dark matter halo is 
1.4xl0 12 M 0 within its virial radius (250 kpc). Warm- 
hot and hot ionized gas dominate the mass accretion over 
large radii with an accretion rate of M=3-5 M 0 yr _1 . 
Nearly 70% of such flow enters the virial radius along 
filamentary structures. This inflowing gas interacts with 
existing halo gas, with its cooling rate and heating rate 
generally balanced over large radii. Only a small fraction 
of the filamentary flows manage to reach regions close to 


3.1. Constructing the Mock Galactic Coordinates 
System in the Simulated Galaxy 

The simulated disk is centered at the point of maxi¬ 
mum baryon density, close to the center of mass of the 
disk. The rotation axis of the simulated disk is aligned 
with the angular momentum vector of the disk cylinder. 
In the mid-plane of the simulated disk, we place a mock 
observer at a distance of Ro=8.0 kpc^Jfrom the disk cen¬ 
ter. This location is used to set up the LSR frame for 
the synthetic observations, which is the reference point 
for the kinematic properties of the sim ulated CGM unless 
otherwise specified. In Section |4.2.3[ we will additionally 
examine whether the observed kinematic properties de¬ 
pend on the location of the LSR by relocating the mock 
observer to seven other locations in the mid-plane at the 
same galactocentric distance. 

Originating from the LSR, we define UVW axes under 
the right-handed convention: U points towards the galac¬ 
tic anti-center and W aligns with the disk rotation axis. 
The direction of a sightline is specified by the galactic 
longitude and latitude (/,6) computed using: 


b 


n ( \ 

— — a,Yccos(w ■ r) 


(i) 


3 We note that Ro=8.0 kpc is slightly shorter than the current 
best estimate of the G alactic distance of the Sun (Ro=8.34±0.16 
kpcjReid et al. |2014 |). We calculate the observable mass fraction 
fm (see Section 14.21 * tor detailed description of / m ) with 8.34 kpc 
and find it only changes by 0.15%, thus the effect of increasing Ro 
by 0.3 kpc is negligible. 


























3 


300K 



0.2 0.4 

f m( V LSR> 0km S_i ) 


0.6 

-l > 


0.8 


0.99 


Fig. 1.— Mass fraction of red-shifted gas (^lsr > 0) in the simu¬ 
lated CGM. The map is in Aitoff projection with galactic center at 
the center of the figure and galactic longitude increasing leftwards. 
Due to the circular motion of the LSR, the sky is heavily red-shifted 
in directions at 0 < l <180 and blue-shifted at -180 < l <0 . 

The diminishing color scales towards higher latitudes indicate the 
existance of both red-shifted and blue-shifted gas. 


and 

^ _ f tv + arccos(tx • r uv ) if 0 < arccos(t; • r) < | 

— 1 7r — arccos(tx • r uv ) if | < arccos(v • r) < 7r, 

( 2 ) 

where r uv is the sightline vector r projected onto the 
disk UV plane, and u, v, w are unit vectors along the 
UVW axes, respectively. 

The mock observer co-rotates with the simulated disk, 
and its circular velocity is approximated by the bulk mo¬ 
tion of baryons within 1.0 kpc of the LSR. To visualize 
the overall gas motions in the simulated CGM with re¬ 
spect to the mock observer, we generate sightlines with 
a spatial interval of (51, 5b) = (1,1) and calculate the 
mass fraction of red-shifted gas (^lsr >0) along each 
sightline. Figure [l] shows the spatial variation of the 
mass ratio in Aitoff projection, with the galactic center 
located at the center of the figure and galactic longitude 
increasing leftwards. Velocity segregation across Z=180 
can be observed: gas at 0 < l <180 is largely red- 

shifted and that at -180 < / <0 blue-shifted. This is 
consistent with observations of the MW’s CGM whose 
velocities are greatly shaped by the Galactic differential 


Wakker & 


rota tion in t he view of local observers (e.g. _ 

Woerden||1991[ ISembach et al.||2003[ |Collins et al. 
2005||Shull et al.||2009|). ~ 


3.2. Kinematic Ranges of the Simulated Disk Rotation 

The strong emission or absorption from the disk gas 
has long plagued observations of the MW’s CGM. In 
most directions the Galactic disk gas appears as low- 
velocity spectral line features that can be efficiently ex¬ 
cluded with certain velocity cutoffs. The adopted veloc¬ 
ity range for the gas associ ated with the MW disk is often 
I^lsrI < 90 — 100 km s" 1 ( Wakker 1991a; Wak ker fc van 
Woerden]|1997 Sembach et al.||2UU3| ~ Lehner et al. 2(J12| 


! 


where rrsR is the velocity with respect to the LSRT 

morejrccurate value is |rr>Ev| < 50 — 75 km s -1 (Wakker 
19 91a[ PPJ12) whe re ^dev is the deviation velocity given 
by’Wakker (1991a) as, 
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Fig. 2.— H I intensity as a function of LSR velocity and galactic 
longitude (l) at b = 0° (top) and b = 70° (bottom). The black 
solid curves represent the T*, = 5 K contour used to set the galactic 
rotation limits (vmin, vmax) in this paper. The green dashed curves 
show the |^dev| = 50 km s -1 dividing lines beyond which the gas 
in the simulated CGM is observable and the horizontal magenta 
lines show the I^lsrI = 100 km s —1 dividing lines. 


^DEV = 


^LSR - V, 
^LSR 


mm 

l^max 


if vlsr < 0 
if v LSR > 0 


(3) 


. The ('LVmn, Vmax ) above is the minimum and maximum 
LSR velocities compatible with the Galactic differential 
rotation in a given direction. This can be obtained either 
by assuming a simple Galactic rotation model or tracking 
the moti ons of neutral hydrogen or molecular clouds i n 
the disk (|Hartmann & Burton||19971 |Dame et al.||2001|) . 

In our synthetic observations, the ( v rnin , v rnax ) is calcu- 
lated from the synthesized H I 21 cm emission line from 
the s imulated disk (following th e description in Draine 
2011| also see |Kim et al.||2014|). We approximate the 
(■ Vrnin: v max ) by the velocities beyond which the H I 
brightness temperature T^ drops below 5 K in an H I 
spectrum. Figure [2] shows two H I position-velocity maps 
at b = 0 and b = 70 , respectively. The 5 K temper¬ 
ature threshold defining the (v m i n , v max ) is marked by 
two black solid curves. As seen in the top panel, the 
(■ VmimVmax ) are very direction-sensitive which indicates 
the strong impact of the circular motion of the LSR on 
the observed radial velocities at low galactic latitudes. 
In the simulated disk, a small warp is found in the +V 
direction bending towards the -W axis and slightly goes 
beyond the surface of the disk cylindrical region (also see 
Figure [7]). This causes the asymmetry of the H I 21cm 
intensities at l = ±90° in top panel. H I intensities be¬ 
low 5 K do not extend significantly, indicating that this 
cutoff encompasses most of the H I 21 cm emission from 
the simulated disk. 

Using (■ Vmin , Vmax), we define two velocity cutoffs to 
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Fig. 3.— The position-velocity(LSR) distribution of the CGM 
(T<10 6 K) when observed from within the simulated disk is shown 
by the shaded region. The mass-weighted LSR velocities of the sim¬ 
ulated CGM gas are shown, with lighter shades indicating a higher 
density of material. The velocity cutoff at |ulsr| — 100 km s -1 
is shown by the horizontal black lines and the simulated CGM 
within this cutoff is lightened. Color-coded symbols show the ve¬ 
locity cent roids of high-velocity gas in the MW ’s CGM from ob¬ 
serva tions (Lehner et al. 2012 Shull et al. 2009 and S embach et al.| 
2003). Open symbols indicate tlie association of the observed high- 
velocity gas with the Magellanic System fM S). 

- Plot generated with the software astroML ( jlvezic et a l. 2014J. 

mask out the low velocity gas in the simulated disk: 
|^dev| > 50 km s _1 , which closely traces the shape 
of the galactic rotation in each direction; and |ulsr| > 
100 km s _1 , which is insensitive to direction but also in¬ 
dependent of disk rotation models. Gas with velocities 
beyond (within) these ranges will be referred as “high 
(low) velocity gas” throughout the paper. In Figure 
[$J the |ulsr| = 100 km s -1 and |udev| = 50 km s -1 
velocity cutoffs are shown as the dotted magenta lines 
and the dashed green curves. The velocity cutoff of 
|ulsr| > 100 km s -1 fails to cover high velocity gas near 
the galactic center which is still associated with the disk, 
while in other directions it is so wide that excess gas in 
the simulated CGM is excluded. 


4. RESULTS 

4.1. Kinematic Similarities: 

Simulated CGM vs. Galaxies 

We first compare the kinematic properties of the sim¬ 
ulated CGM (T<10 6 K) and that observed for the MW. 
For each sight line through the simulated CGM, we cal¬ 
culate the mass-weighted LSR velocity (ulsr) over all of 
the simulated gas cells as shown in Figure [3| Observa¬ 
tions of high velocity gas in the MW’s CGM from various 
metal lines are shown as color-co ded symb ols (e.g., Si_II, 



Magellanic System is excluded (open symbols at gener¬ 
ally I^lsrI >200 km s _1 ), there is overall consistency in 
the gas kinematics between the simulated CGM and the 
MW’s CGM. The velocity spread of the gas in the simu¬ 
lated CGM is narrower than that of the actual MW. This 
is most likely due to the mass-weighted algorithm which 
smooths out high-velocity low-density features. Though 
the simulated MW does have small satellites, they are 
not at the scale of the Magellanic System and contribute 


Fig. 4.— The velocities of the CGM (T<10 6 K) with respect 
to the galaxy systemic velocity as a function of impact parameter. 
The gray shades show the line-of-sight mass-weighted velocities of 
the gas in the simulated CGM when observed from three external 
views: edge-on, 45 and face-on. The three horizontal dash lines 
indicate the median values of the gas velocities from the edge-on 
view (bottom), 45 (middle) and face-on view (top). Color-coded 
symbols show the gas centroid velocities o f different transitio n lines 
from the COS-Halos survey (|Tumlinson et al.||2013| |Werk et al.| 
|2013|). Multiple data points at the same impact parameters rep- 
resent the multiple component structures along the corresponding 
line of sight which results from individual Voigt profile component 

analysis. __ 

- Plot generated with the software astroML ( |Ivezic etTa l. 2014). 


negligible mass flux (Joung et al. 2012). In addition, 
three of four satellites in the simulation at z^0 are at ve¬ 
locities <200 km s -1 with respect to the galaxy systemic 
velocity. 

Synthetic observations of the simulated CGM from ex¬ 
ternal views are also conducted so that we can make 
connections with other L~L* galaxies. We observed the 
simulated CGM (T<10 6 K) with three disk inclinations: 
0 (face-on), 45 , and 90 (edge-on). From each inclina¬ 
tion, we calculate the line of sight velocities for 500 x 500 
sightlines on a square grid that spans -250 kpc to 250 
kpc (with intervals of 1 kpc) from the galactic center. 
In Figure [4j the line of sight velocity Sv with respect to 
the galaxy systemic velocity as a function of impact pa¬ 
rameter is plotted as the gray shades for the simulated 
CGM. The observation symbols are from the COS-H alos 
survey (Tumlinson et ai.l 2013 Werk et al. 2013) and 
are similarly color and shape-coded as Figure [3] show¬ 
ing the observed velocity centroids of different species 
(i.e., H I, Mg II, C II, C III, Si III and O VI) in the 
CGM of L^L* galaxies at z^0.2. The velocity disper¬ 
sion of the simulated CGM is a ~ 53 km s -1 , which 
is smaller than that derived from the C OS- H alos survey 
(a ~ 85 km s -1 from the metal lines by Werk et al. 2013 
and a ~ 200 km s -1 in the H I by Tumlinson et al 


2013). Some of the difference could be the lack of AGN 


feedback in the simulation, but more important effect is 
likely that only one simulated galaxy is examined in our 
synthetic observations while in the COS-Halos’ sample 
galaxies may slightly vary in mass, gas content and so 
on. High velocity outlie rs can be seen from the H I data 
(Tumlinson et al. 2013), that may come from potential 
fast moving satellites m the CGM of the L~L* galaxies 
or non-relevant cold gas parcels beyond the CGM but 
intersected by the sightlines. Overall, the kinematics of 



























































5 




Fig. 5. — Observable mass fraction / m of the simulated CGM 
with velocity cutoffs I^devI > 50 km s -1 (top) and I^lsrI > 
100 km s —1 (bottom). Figures are in the same configuration as 
Figure [l] Sightlines that contain primarily high-velocity CGM are 
the darkest, while regions where the CGM is largely obscured are 
the lightest. 


gas in the simulated CGM when viewed externally is in 
good agreement with the observations of the CGM of 
other L L* galaxies. 

When viewed from different inclinations, the mass- 
weighted velocity of the CGM gas shows asymmetric dis¬ 
tributions. This is indicated by the dashed lines in Figure 
[4| which are the median values of the velocity distribu¬ 
tions for edge-on (bottom), 45 degree (middle) and face- 
on (top) views. We find a bulk motion of the simulated 
CGM with respect to the entire galaxy in a direction par¬ 
allel to the disk plane with Sv of 30 km s -1 . This bulk 
flow therefore contributes to the line of sight velocities of 
the simulated CGM from the edge-on view, while obser¬ 
vations with th e face-on view wi l l large ly miss it. This is 
consistent with Kacprzak et al. (2010) who showed that 
the Mg II absorption line profiles m their sample favor 
one side of the galaxy systemic velocity, suggesting that 
the observed velocity spreads may depend on galaxy in¬ 
clinations. 


4.2. Gas Observability in the Simulated CGM 
4.2.1. Overall Gas Observability 

We first calculate the mass fraction and spatial distri¬ 
bution of high velocity gas in the simulated CGM. For a 
given sightline (/, 5), the mass fraction / m of high veloc¬ 
ity gas along a path length of 250 kpc is computed. The 
top (bottom) panel in Figure [5] shows the variation of fm 
across the simulated galactic sky with a velocity cutoff 
of |^dev| > 50 km s -1 (|^lsr| > 100 km s _1 ). Close 


TABLE 1 
Observability 


Components 

(i) 

Norm. Mass 
(2) 

fm,v dev 

(3) 

fm,v lsr. 

(4) 

CGM (r<250 kpc) 

1.00 

0.54 

0.45 

|b| >10“ 

- 

0.44 

0.40 

6 >20° 

- 

0.33 

0.35 

T < 10 5 K 

0.08 

0.62 

0.55 

10 5 < T < 10 6 K 

0.70 

0.54 

0.44 

T > 10 6 K 

0.22 

0.53 

0.44 


Note. — Column 2: gas mass normalized to the total 
mass of the halo; column 3: mass fraction of high velocity 
gas with I^devI > 50 km s — 1 ; column 4: gas mass frac¬ 
tion with I^lsrI > 100 km s -1 . Row 2: overall gas ob¬ 
servability in the simulated CGM; row 3: observability of 
halo gas at |6| > 10 and \b\ > 20 ; row 4: gas observabil¬ 
ity within temperature ranges T < 10 5 K, 10 5 < T < 10 6 
K and T > 10° K, respectively. 


to Z±90°, large patches with f m -A 1 are observed. This 
indicates that most gas along corresponding sightlines 
move at high velocities which are o bserv able to the mock 
observer. As mentioned in Section [T2j the velocity cut¬ 
off I^lsrI > 100 km s -1 is generally too broad for most 
directions, resulting in more obscuration of the simulated 
CGM. This reflects in Figure [5] (and Table [l]), where f m 
with |^ L sr| >100 km s -1 is generally lower than that 
with I^devI > 50 km s -1 . At the bottom of both panels 
exists a large orange patch, resulting from high velocity 


SN- driven outflows from the simulated disk (Joung et al. 


2012). This outflowing gas can also be found in the phase 


diagrams (Figure |6| at temperature above 10 6 K. 

Integrating f m over the whole spherical volume 
(r <250 kpc), we obtain the total mass fraction of high 
velocity gas for different CGM components defined b; 
temperature ranges and galactic latitudes. In Table 
column 2 shows the mass of each component norma’ 
ized to the total gas mass of the simulated CGM. The 
/m,v DEV (/m,v L sr) i n column 3 (4) indicates the mass 
jfraction of high velocity gas with the velocity cutoff 
|^dev| > 50 km s -1 (|^lsr| > 100 km s _1 ). We find 
that generally only ~0.5 of the CGM mass is carried by 
high velocity gas; this means that nearly half of the mass 
of the simulated CGM is obscured from observations due 
to velocity cutoffs. This is one of the key results of this 
paper. 


Specifically, the mass fractions at |6| > 10° and \b\ > 
20 are computed since low latitude regions in the actual 
MW’s CGM are usually avoided by absorption line exper¬ 
iments due to Galactic emission or absorption. We find 
that the / m?VDEV (/ m ,v LSR ) decreases to ^0.4 at |6| > 10 
(~0.3 at \b\ > 20 ), indicating that the observations of 
the actual MW’s CGM may miss more than half of the 
gas mass. It is notable from Figure [5] that the simulated 
southern sky is covered by more high velocity gas than 
the northern sky. The MW (and other actual galaxies) 
may have this type of imbalance but it is biased by the 
presence of the Magellanic System in the southern sky. 


4.2.2. Gas Observability at Various Phases and 
Galactocentric Radii 
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Fig. 6.— Phase diagrams of the simulated CGM in galactocentric shells at [20,30], [50,60], [110,120], [200,210] kpc. High velocity gas 
moving at |udev| > 50 km s -1 (i.e., observable) is shown by horizontal red strips, while that with |udev| < 50 km s -1 is marked by 
vertical green. 


In order to study if velocity cutoffs preferentially select 
gas at certain phases, we calculate the mass fraction of 
high velocity gas for gas with temperature T < 10 5 K, 
10 5 K < T < 10 6 K and T > 10 6 K. The last row in 
Table [T| shows that the mass fraction of high velocity gas 
remains nearly 50% for gas at different phases in spite of 
their distinct contributions to the total CGM mass. We 
compare the phase properties of the high velocity and low 
velocity gas in the simulated CGM in Figure |6j and find 
that they occupy similar regions in the phase space (see 
below). This indicates that observations of the MW’s 
CGM using various metal lines from different ionization 
states may miss similar amounts of low velocity gas no 
matter which temperature ranges they probe. 

The H I HVCs in the actual MW’s CGM are generally 
observed within 0.06 R v ir,MW while the COS-Halos team 
finds CGM out to 0.55 R v ir,cos Q To study if velocity 
cutoffs preferentially select high velocity gas at certain 
radii, we decompose the volume of the simulated CGM 
into a set of galactocentric shells, each with a thickness 
of Ar=10 kpc. We separately calculate the mass frac¬ 
tion of high velocity gas within each shell, and find that 
f m fluctuates very little from r=10 kpc to r=250 kpc, 
with an average of 0.54T0.03 for |^dev| > 60 km s -1 
(0.44T0.03 for |^lsr| > 100 km s _1 ). Such constant f m 
throughout the simulated CGM can be explained if the 
gas velocities of the simulated CGM are re lativ ely contin¬ 
uous along radial directions (see Section 4.4 for further 
discussion). 

Figure [6] shows the phase of high and low velocity gas 
within the different galactocentric shells. Only the re¬ 
sult with I^devI > 60 km s -1 is shown since that with 
|^lsr| > 100 km s -1 is similar. Though the high and low 
velocity gas in the simulated CGM do generally occupy 
similar phase space, the high velocity gas is more ex¬ 
tended to the high-temperature, low-density region (the 
top left corner) which correspond s to hig h ve locity, hot 
outflows in the simulation (Joung et al. 2012, PPJ12). 
There is also some high density, low temperature gas in 
the shells between 50-120 kpc that is only found in the 
high velocity regime. This corr espo n ds to gas st ripped 
from recently accreted satellites ( Joung et al .|2012 ) and is 


4 We adopt R v ir,MW as 250 kpc, same as the value we used for 
our synthetic obs ervations; and R v i r .cos ~300 kpc is given by the 
COS-halo teams (|Werk et al.|2014| >. 


a negligible fraction of the total CGM mass. We also ex¬ 
amined the simulated CGM at different metallicity bins 
ranging from 0-0.1 Z© to >1 Z©, and find that the mass 
fractions of high velocity gas are similar (/ m ~0.6) in 
spite of the distinct mass distribution at each bin. These 
results suggest that we are not biased towards high veloc¬ 
ity gas (T < 10 6 K) having specific properties that may 
represent distinct origins. 

Within each galactocentric shell, high velocity gas can 
generally be found at directions towards /±90 which 
is similar to the spatial distribution of the overall case 
shown in Figure [5] This is to say, the velocity fields 
in the simulated CGM are closely correlated in space 
throughout the volume. This spatial coherence is in large 
part due to the circular motion of the LSR, which shifts 
halo gas around /±90 to high velocity by adding extra 
components to the gas radial velocity with respect to 
the LSR. On the other hand, such large scale velocity 
correlation is also aided by the presence of large scale 
gas flows, such as feeding filaments and polar outflows. 


4.2.3. Varying the Location of the Mock Observer 

Gas observability in the simulated CGM discussed 
above can be biased since the mock observer is so far 
fixed at one location in the disk mid-plane. Local activi¬ 
ties such as turbulence and SN-driven outflows may cause 
non-negligible effects on the corresponding synthetic ob¬ 
servations. In this section, we vary the location of the 
mock observer in the disk mid-plane and study if the 
change of the reference LSR affects the observed gas kine¬ 
matics in the simulated CGM. 

We relocate the mock observer to seven other places in 
the simulated disk with an angular separation of 0 = 45 , 
where <f> is the polar angle on the fiducial UV plane de¬ 
fined in Section [3] with +U being the polar axis. Each 
location is 8.0 kpc from the galactic center and the ve¬ 
locity of a given LSR is again averaged over the baryons 
in its local neighbourhood in order to closely represent 
the circular motion of the disk rotation at that location. 


We follow the same procedures described in Section 4.2.1 
and calculate the mass fraction of high velocity gas in the 
simulated CGM from corresponding new locations. We 
only apply |^lsr| > 100 km s -1 in our calculations but 
note that the values with |^dev| > 60 km s" 1 will be 
similar considering the results in Section 4.2. 1| 

Table [2] shows the mass fraction fm,v LSK of high veloc- 
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TABLE 2 

Variation of Observability 











0 

45 

90 

135 

180 

225 

270 

315 

fm,v lsr 

0.45 

0.39 

0.55 

0.70 

0.67 

0.64 

0.53 

0.45 


Note. — Mass fraction of high velocity gas as a function of the 
locations of the mock observer. Velocity cutoff is set as |ulsr| > 
100 km s -1 . Each location is 8.0 kpc from the galactic center, and 
</> = 45 apart where 4> is the polar angle on the UV plane defined 
in Section[3] 


ity gas in the simulated CGM as observed from new lo¬ 
cations. fm, v lsr = 0-45 at 0=0 corresponds to the value 
from the original LSR, which has been shown in Table 
[l] The gas observability in the simulated CGM varies 
Rom 0.39 to 0.75 (with a mean of 0.55) as the angular 
locatio n of the mock observer changes. As indicated in 
Section 


4.1 


a net bulk flow (Sv ^30 km s -1 ) of the sim¬ 
ulated CGlVl with respect to the galaxy is observed in 
our simulation. This may largely cause the variation of 
fm,v lsr si nce the calculations of the velocities of the LSR 
only consider the bulk motions of local neighbourhood. 

We examine the spatial distribution of high velocity 
gas in the simulated CGM as observed from these eight 
different locations, and find similar high velocity patterns 
as have been seen in Figure [5] Except for the hot out¬ 
flowing gas towards the poles, most high velocity gas is 
observed in directions around /±90 , indicating that gas 
velocities at low galactic latitudes are greatly shaped by 
the circular motion of the LSR in the disk (see also Figure 
[T]) . In all cases, we still find that nearly half of the CGM 
mass is carried by low velocity gas, which is obscured due 
to velocity cutoffs in spite of different locations of mock 
observers in the disk plane. 

4.3. H I Observability in the Simulated CGM 

In this section, we calculate the distribution and cov¬ 
ering fraction of high velocity H I in the simulated CGM. 
We only present the H I behaviors with velocity cutoff 
|^dev| > 50 km s -1 , with the knowledge that the results 
with |^lsr| > 100 km s -1 are similar according to Figure 
[5j For a given sightline, we compute the column density 
of high velocity H I (Nhi) by summing up the H I retained 
by all the simulated gas cells with |^dev| > 50 km s -1 
along the corresponding path (^250 kpc). We also use 
different path lengths ranging from 50 to 300 kpc and 
find little variations among those Nhi values we have ob¬ 
tained. This is because most of the H I resides at radii 
close to the disk except for satellites, an increase in path 
length will not necessarily gather more H I. 

Figure [7] shows the spatial distribution of high veloc¬ 
ity H I in the simulated CGM. The high Nhi clumps 
near (/, b) = (—90 , —30 ) are the warp of the disk bend¬ 
ing downwards which was not in clud ed in the cylindrical 
disk region as defined in Section 3.2 (see also |Fbrnandez| 
et al.|2012 ). This uneven distribution of H I at / = ±90 
causes the asymmetric H I 21 cm emission in Figure [2] 
At higher elevation above the disk plane, elongated struc¬ 
tures can be seen which may relate to gas-rich satellites 
or accreted gas along filamentary streams in the simu - 
lated CGM (Joung et al.||2012| |Fernandez et al.l|2012|). 
Besides these distinct dense structures, most of the sim¬ 



FlG. 7.— Column density of high velocity H I (|udev| > 
50 km s -1 ) in Aitoff projection. The galactic center is located 
at the center of the figure, and galactic longitude l increases from 
left to right. 


ulated galactic sky is covered by very diffuse H I. 

Figure [8] shows the covering fractions of high velocity 
H I and the H I without velocity cutoffs in the simu¬ 
lated CGM at different column sensitivities. Regions at 
\b\ < 20 are purposely omitted in order to link to the 
actual observations of the MW’s CGM. A comparison be¬ 
tween the black solid line (without velocity cutoffs) and 
the green-dash lines (with velocity cutoffs) indicates that 
the covering fraction of the high velocity H I is reduced 
to only 0.29 (from 0.48) at Nhi > 10 18 cm -2 . More 
generally, the covering fractions of the high velocity H I 
are ^0.2 less than those without velocity cutoffs at a 
given column density sensitivity. Color-coded symbols at 
Nhi ^ 10 18 cm -2 mark the measured covering fractions 
of HV Cs in the a ctual MW’s CGM: 0.18 at Nm > 2xl0 18 
cm" 2 jWakkerl 199 lap 0.37 at Nhi > 7 x 10 17 cm -2 
( Murphy et al.||19 95), and 0.37 at Nhi > 8 x 10 17 cm -2 
(Lockman et al. "2002). The consistency between the 
measured HI HVCs and the simulated high velocity H I 
column density indicates that a considerable fraction of 
low velocity H I m ay r emain hidden in the actual MW’s 
CGM (see Section [543] for a detailed discussion of poten¬ 
tial missing low velocity H I complexes). 


4.4. Additional Mass Obscuration in the Outer Halo 

Apart from the mass obscuration due to foreground low 
velocity disk gas, high velocity gas in the outer halo may 
be additionally blocked by nearby halo gas if it moves at 
similar velocities. In another words, even at high veloc¬ 
ity regimes spectral lines from the gas in the outer halo 
may blend with signals from nearby gas if line separa¬ 
tions (i.e., velocity differences) are too small to detect. 
Recognition of these close line profiles will depend on the 
resolution of observations. 

To quantify this issue, we extract two data sets from 
the simulation which represent an inner halo (r <30 kpc) 
and an outer halo (50< r <150 kpc), respectively. Maps 
of high velocity gas in both halos in Aitoff projection 
show similar patterns as those in Figure [5] For a given 
sightline, we calculate the velocity differences between 


5 No latitude limitation |6| > 20° 
(1991a) . 


was adopted by 


Wakker 

















































Fig. 8. — H I covering fractions vs. H I column sensitivities in 
the simulated CGM (|6| > 20°). Black solid line: H I without any 
velocity cutoffs; green dashed line: H I with |r>DEv| > 50 km s -1 . 
Color-coded symbols: observed H I HVCs covering fractions in the 
actual MW’s CGM b vlWakker ( 1991a]) , |Lockman et al.| (|2002|) and 
|Murphy et al.||l995| ). 


the high velocity gas in the inner halo and that in the 
outer halo. This is separately performed for red-shifted 
and blue-shifted gas. Gas in the outer halo is observable 
only if: (1) its velocity is high enough to survive the 
J^devI > 50 km s -1 or |^lsr| > 100 km s -1 cutoffs; (2) 
its velocity differs from that of the gas in the inner halo 
by a number Av which is the chosen minimum detectable 
line separation. 

Table [3] shows the observable mass fraction of high ve¬ 
locity gas in the outer halo that satisfies both criteria 
under a set of minimum detectable line separations Av 
ranging from 0 km s -1 to 30 km s -1 . If the spectral res¬ 
olution is high enough (0< Av <10 km s -1 ) so that gas 
at similar velocities in the inner and outer halo would be 
distinguishable, we find that the observable mass frac¬ 
tion in the outer halo is 0.47 at |vdev| > 50 km s -1 
(0.43 at I^lsrI > 100 km s _1 ). However, the observable 
mass fraction in the outer halo decreases as the spectral 
resolution falls off (i.e., Av becomes larger). For exam¬ 
ple, it is significantly reduced by a factor of two (from 
0.47 to 0.23) if Av = 30 km s -1 . This indicates that less 
than 50% of gas is accessible to local observers due to 
a combination effect of velocity cutoffs and obscuration 
by nearby halo gas. In particular, this has a potential 
influence on the mass estimate of the ionized HVCs by 
Lehner et al. (2012) who compared the covering fractions 
from a stellar sample with that from a QSO sample, and 
constrain the clouds to be within 5-15 kpc. 


5. DISCUSSION 

The CGM is key to our understanding of galactic ac¬ 
cretion and feedback processes. It has been claimed that 
the CGM can solve the “missing baryon problem” Jrom 
studies with both obs ervations (e.g., |Werk et al.||2014|) 
and simulati ons (e.g., Cen & Ostriker " 2006 ” Sommer 


Larsen 2006). In this section, we re-evaluate the total 
mass of the MW’s CGM taking into account the ob¬ 
scured mass due to the velocity cutoffs, and compare 
the obscuration-corrected MW’s CGM with the CGM of 
other L^L* galaxies. In the end of this section, we dis¬ 
cuss potential observations of obscured low velocity gas 


TABLE 3 

Gas observability in the outer halo 
(50 < r < 150 kpc) 


Min. Sep. 

0 

10 

20 

30 

Av 

km s -1 

km s -1 

km s -1 

km s -1 

fm,v dev 

0.47 

0.39 

0.30 

0.23 

fm,v lsr 

0.43 

0.31 

0.23 

0.16 


Note. — Av means a minimum detectable line 
separation in spectroscopic observations; /m,t; DEV and 
fm,v LSR are the observable mass fractions of gas in 
the outer halo at |udev| > 50 km s -1 (|ulsr| > 
100 km s -1 ) with the corresponding Av. 


in the MW’s CGM through non-kinematic methods. 

5.1. Obscuration-corrected Mass of the MW’s CGM 

Cold gas (T< 10 4 K) in the actual MW’s CGM has 
been broadly studied with H I 21 cm observations, while 
the warm (^10 4-5 K) and warm-hot gas (^10 5-6 K) 
has been revealed by UV absorption lines in the spectra 
of background sources (stars or QSOs) and Ha emission. 
Mass estimates of this gas focus on available high velocity 
lines in spectra, with little evaluation of the mass contri¬ 
bution from low velocity obscured halo gas. Thus in this 
section we apply the observable mass fraction f m = 0.5 to 
current mass estimates of the MW’s CGM and obtain 
an obscuration-corrected mass for CGM gas at various 
phases. We do not apply the correction to gas from the 
Magellanic System as this system is moving at hig h ve¬ 
locity relative to the MW (Kallivayalil et al. 2013) and 
the vast majority of the gas is likely to be observable 
based on prediction s from the models (Be sla et al.||2010 


Connors et al.|2006 ). This type of large accretion event is 


also not presently occurring in our simulated galaxy halo, 
and i s rare in observations of MW-like galaxies (Tollerud 
et al. 12011 ). 


Cold gas in the MW’s CGM is usually in the form 
of discrete H I structures called HVCs. Galactic HVCs 
have a total mass of Mhi ^3x 10 7 M q , excluding the 
HVCs associated with the Magellanic System. This in¬ 
cludes 13 H I complexes within 10 kpc of the disk and 
the remaining small HVCs placed at a general distance 
of 10 kpc (PPJ12). Considering the 50% mass obscu¬ 
ration, the total mass of Galactic HVCs is therefore 
corrected to M^[ r ^6xl0 7 M©. This corrected H I 
mass in the MW’s CGM without the Magellanic Sys¬ 
tem is consiste nt with observations of M31’s halo gas 
(Thilker et al. 2004). If we include the mass from the 
high velocity gas of the Magellanic System (excluding 
the Magellanic Clouds themselves),the mass inc reases 
by Mht = 3 x IQ 8 (d/55 kpc) 2 M© (|Putman et al. 2003 
Bruns et al.||2005| ). 


There are recent mass estimates of warm and warm- 


hot gas in the MW’s CGM. Lehner & Howk (2011) found 
that the ionized HVCs reside within 5-15 kpc using the 
constraint from halo stars. The total mass of warm and 
warm-ho t gas is thus M wa , rrn ^l.lxl0 8 (d/12 kpc) 2 M© 
(see also Shull et al.j 2009). If the observable mass frac¬ 
tion / m =U3" is applied, the mass of warm and warm-hot 
gas will be M^ a r r r m ^2.2 xlO 8 (d/12 kpc) 2 M©. This mass 
is almost certainly a lower limit given the consistent find¬ 
ing of warm and warm-hot gas at large radii in the CGM 
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of other galaxies (e.g., Tumlinson et al.|2011 Chen et al. 


2010 Werk et al. 2014|). Some of the lVlVV's COM at large 
radii may b(Ttoo diffuse to detect from within the disk 
using bac kgro und QSOs. One should also consider that 
in Section [474] we show that distant gas in the CGM can 
be obscured by nearby gas if they are moving at similar 
velocities. This shadowing effects from high velocity gas 
in the inner halo may cause an additional 50% of mass 
obscuration in the outer halo in observations with a spec¬ 
tral resolution of >30 km s -1 . Thus the mass estimate 


by Lehner et al. (2012) shall be interpreted with caution 
as they do not take into account the local halo shadow¬ 
ing effects. We do not attempt to correct for this addi¬ 
tional obscuration here or calculate a mass for gas that 
does not have a distance constraint. A recent estimate 
of the total mass in warm and warm-hot gas from the 
Magellanic System (excluding the M agellanic C louds) is 
M wa , rTT 1 .MS ~l-5xl0 9 (d/55 kp c) 2 M 0 ( |Fox et ar||2014 ). 


Miller & Bregman| (2014) estimate the mass ol hot 
gaT[T>To t) K) in the MW s CGM to be 4.3 xlO 10 M 0 
within 250 kpc (or 3.8xlO 9 M 0 within 50 kpc) using 
O V II an d O VIII li nes through X-ray observations (see 
also [Gupta et al.|[2012]). This mas s is consistent w ith 
no n-spe ctroscopic observations (see [Anderson & Breg- 
man|2010| for a detailed discussion). Given the low veloc¬ 
ity resolution of X-ray observations, the mass estimates 
of hot halo gas are usually based on line profile fitting to 
exclude local emission components rather than applying 
velocity cutoffs. Thus our mass obscuration correction is 
not applicable to their mass estimates. In Figure [6| we 
also show that the gas that would be at highest veloci¬ 
ties and potentially directly detectable is the most diffuse 
and hot and the most difficult to detect. 

5.2. Is the MW’s CGM an Outlier? 

In this section we compare the MW’s CGM to that of 
other L^L* galaxies for which large amounts of warm 
and warm-hot (10 4-6 K) have been found. We start by 
looking at the mass estimates and then, given the model¬ 
ing required to obtain total masses, assess the difference 
in terms of the column density. As shown in Figure [4, ex¬ 
ternal galaxies do not have the same obscuration effect; 
the observations of background QSOs often show gas in 
CGM very close to the systemic velocity of the galaxy. 

5.2.1. CGM Mass 


Stocke et al. (2013) found ~10 10 M 0 of CGM (T-10 4 
K) within the virial radius of late-type galaxies at z< 
0.2 (with 75% of the mass contai ned in the inner half 
of the virial radius). jjWerk et al. (2014) also suggest a 
conservative estimate of 2xlO iU M 0 for the 10 4 K CGM 
gas within 160 kpc of the L^L* galaxies at z^0.2. This is 
^3xl0 9 M 0 within 50 kpc of the disk and ^4xl0 8 M 0 
within 15 kpc based on the gas surface density profile 
they ad opt. The value is comparable to the MW estimate 
(Section [5T]) within small radii, but the MW value is still 
low and the external galaxy results indicate there is most 
likely a large amount of gas undetected at large radii. 

Since the COS results are for galaxies at somewhat 
higher redshift, we examined the simulated CGM at dif¬ 
ferent redshifts to assess if any discrepancy could be due 
to gaseous halo evolution with time. Negligible varia¬ 
tion is found for the mass of the simulated CGM from 


z=0.3 to z=0.0 (within a_fixed radius of 250 kpc). This 
is consistent with Stinson et al. (20l2| who examined the 
CGM gas content m an L=1”C* star-forming galaxy, and 
found little variation of CGM mass within the virial ra¬ 
dius since z=0.5. Thus, the difference in redshift of the 
external galaxies and the MW is not expected to greatly 
affect the CGM mass. 

There are some remaining issues with the CGM mass 
estimates for external galaxies. The warm CGM pres¬ 
sure that Stocke et al. (2013) derive (P/k~10 cm -3 K) 
is an order of magnitude lower than that typically found 
from observations of clouds in the halo of the Milky Way 
(Fox et al. 2010 Hsu et al. 2011). This may be due to 
uncertainties related to the photoionization modeling. In 
addition, onejff the mass estimates for the 10 4 K CGM 
by Werk et al. (2014|) is based on the assumption that 
the nni along each sight line is representative of all of the 
gas at a given impact p arameter. However, observations 
of the MW’s CGM (e.g., Sembach et al.|2003l S hull e t al. 
2009) and external galaxies (e.g., I Werk et al.||2014|) in- 
dicate that multiple discrete absorbers exist along sight¬ 
lines. Our synthetic observations from external views 
and from within the simulated galaxy also suggest mul¬ 
tiple absorbe rs at disc rete velocit ies. Although the other 
method that Werk et al. (2014|) applied takes into ac¬ 
count the volume tilling facto”this mass estimate comes 
with large uncertainty since it is largely dependent on 
the calculation of Nhi and tiri- 

Thus far, we cannot confirm that the MW is an out¬ 
lier in CGM mass. Other potential factors, such as the 
completeness of observations and the gas density profile 
of the CGM, need to be more fully considered. 

5.2.2. CGM Column Density 

Column densities of metal lines detected in the MW’s 
CGM are general ly lower than tho s e of e xternal galax¬ 
ies. For example, Tumlinson et al. (2011) found a typi¬ 
cal value of log Nqvi=14.5 for 27 out of 30 star-forming 
galaxies with sSFR>10 -11 yr -1 , which is higher than 
the mean value of log (Novi)=14.0 (at sSFR^2-6x10 -11 


yr x ) for the MW’s CGM (Sembach et al. 2003). To 


evaluate this discrepancy in column density, variations 
of path length from observations of the MW’s CGM to 
those of external galaxies need to be considered as well 
as the velocity obscuration effect. 

As an approximation, we compute the column densities 
of warm-hot gas N^h in the simulated halo. We calculate 
the ratio of N^ h /N^ h , where “ex” indicates column den¬ 
sities of the warm-hot CGM as observed from external 
views (edge-on, 45 and face-on), and “in” for that with 
an internal view as described in Section [ 3 J For N^ h , we 
sum up all the material along the sighthnes and do not 
separate high velocity and low velocity halo gas so that 
the influence of path length difference independent of 
disk obscuration effects can be evaluated. We find that 
the ratio of external to internal column density varies 
from 0.9 to 5.0 with a mean of 2.34. This implies that 
the column densities of CGM in external galaxies are in¬ 
trinsically higher by a factor of ^2.0 than those of the 
MW because of the difference in path length. Taking 
this into account, we can correct the column densities of 
the MW’s CGM to log (Novi)=14.30, similar to other 
L^L* galaxies. This is a lower limit given that the low 
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Vlsr (km s' 1 ) 

Fig. 9. — Left: the WA complex from the HVC catalog of (Wakker & van Woerden (19911) is shown in ellipses. The size of the ellipse 
indicates the sky coverage of each WA HVC, but does n ot correspond to the shape of the observed WA HVCs on the sky. The L2+L5 
(Peek et al.|'2009) and +WQ3 LVHCs (|Saul et al. 2012"} are indicated by squares and triangles respectively, although their sizes are not 
accurately indicated. The color of the symbols indicates the LSR velocities. Right: the mass-weighted LSR velocities of the simulated 
CGM of the same patch of sky. The velocities shown in both panels are in the same color scale as is indicated by the color bar. A velocity 
gradient can be clearly seen both in the MW and the simulated CGM. 


velocity gas obscured by MW’s disk could make an addi¬ 
tional contribution to the average column density. Thus 
we find the MW’s halo detected in O VI is not an outlier 
in column density. 

5.3. Obscured Low Velocity H I complexes 

In this section we examine the obscured halo gas in 
the context of complexes of HI gas that may exist in the 
Milky Way halo. These complexes would be low velocity 
halo clouds (LVHCs), or halo clouds that have LSR ve¬ 
locities indistinguishable fro m those in the MW disk, but 
are o therwise akin to HVCs (Wakker 1991b| Peek et al. 
2008). These LVHCs may represent bulk nows within 
the MW halo that have been missed by standard HVC 
studies due to our location in the disk. 

There is some evidence for MW LVHCs from far- 
infrared (FIR) studies of low velocity clouds that find 
a small n umber have undetectable FIR emi ssion, similar 
to HVCs (JPeek et al.||2009| |Saul et al.||2014|). This lack 
of FIR emission in HVCs is likely due to some combina¬ 
tion of low dust content and weak dust h eatin g from the 
interstellar radiation field. Peek et al. (2009) examined 
the ratio of dust FIR emission to H I gas column den¬ 
sity and found several potential L VHC s, in particular L2 
and L5 noted in Figure [9| Saul et al. (2014) examined 
this ratio for a larger population of low velocity clouds 
and identified an entire complex of potential LHVCs in 
a region o f sky near a HVC complex known as the WA 


complex. Saul et al. ( |2012| ) referred to this potential 


LVHC population as the~fWQ3 clouds. As can be seen 
in the left panel of Figure [9| there is a position-velocity 
link between the +WQ3 clouds, the WA complex, and 
L2 and L5, indicating these clouds represent a bulk flow 
that was not previously recognized in halo gas studies. 

To examine if we see similar LVHC features in the sim¬ 
ulated CGM, we plot the LSR velocities of gas in the 
simulated CGM for the same patch of sky as the clouds 
discussed above in the right panel of Figure [9j We find a 
velocity gradient in the simulated galactic sky, which is 
consistent with the observed pattern from complex WA 


to the L2, L5 and +WQ3 LVHCs. This implies that 
those regions with low ulsr towards the galactic center 
and anti-galactic center in the simulated halo have pop¬ 
ulations of LVHCs, similar to what the MW FIR studies 
seem to be finding. Since the low ulsr regions (i.e., the 
low f m regions in Figure J5| are spatially coherent, we 
expect there are additionalTarge LVHC complexes in the 
MW’s CGM which have been missed entirely due to the 
velocity cutoffs provided by the disk emission. 

6. CONCLUSIONS 

When we observe the MW’s CGM, the low velocity 
gas is missing and we are biased towards high velocity 
gas due to disk contamination. To assess the amount of 
obscured low velocity halo gas and its contribution to the 
MW CGM’s baryon mass, we conduct synthetic observa¬ 
tions of a simulated MW-mass galaxy. We embed mock 
observers inside the simulated disk and observe outward 
to examine gas kinematics as shaped by galactic rotation. 
We also observe the simulated CGM externally for three 
different disk inclination angles. The synthetic observa¬ 
tions of the simulated CGM with inside-out views and ex¬ 
ternal views show comparable kinematic properties with 
the MW’s CGM and that of other L ~ L* galaxies. 

Our main results can be summarized as follows: 


1. The mean mass fraction of high velocity (observ¬ 
able) gas in the simulated CGM is 0.55 (Figure [5] 
Table lj as is seen by the mock observer in the sim- 
ulateadisk. This value fluctuates from 0.39 to 0.75 
as the location of the mock observer changes. 


2. Large spatially-coherent velocity fields are seen in 
the simulated CGM, where most of the high veloc¬ 
ity gas is at low galactic latitudes and the low veloc¬ 
ity gas generally exists at higher galactic latitudes. 
This suggests that large H I complexes moving at 
low velocities in the MW’s CGM may be missed by 
observers located at solar radius. 
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3. Given the similarity between the simulation and 
the MW’s CGM kinematics, we find that the 
obscuration-corrected mass of the MW’s CGM (T< 
10 6 K) is 2.8x10 s Mq excluding the Magellanic 
System. This gives the lower bound of the CGM 
mass since current mass estimates are limited to 
the inner halo where reliable distance constraints 
can be found. Although the corrected mass of 
the MW’s CGM is below the CGM mass of other 
L ^ L* galaxies, we find similar column densities 
for O VI lines after correcting for the difference 
in path length between observations of the MW’s 
CGM and that of external galaxies. 

4. Obscured low velocity gas and observable high ve¬ 
locity gas in the simulated CGM closely overlap in 
phase space, indicating that this 50% obscuration 
in mass affects observations of gas at most phases 
at T < 10 6 K (Figur e jgb . This is to say, observa¬ 
tions of the MW’s CGM using various metal lines 
at different ionization states generally miss 50% of 
CGM gas because it moves at low velocities. In 
addition, no significant difference in metallicity is 
found between the low velocity and high velocity 


gas. 

5. Breaking the simulated CGM into a series of galac- 
tocentric shells, we find that the observable mass 
fraction does not vary with galactocentric distance 
(Figure [6]). The spatially coherent distribution of 
high velocity gas in the inner and outer halo may 
significantly reduce the gas observability in the out¬ 
skirts of the MW’s CGM. This happens if distant 
high velocity halo gas moves at velocities close to 
that of the nearby halo gas so that their line sepa¬ 
rations in spectra are too small to recognize. 

We thank Greg Bryan for useful discussions on the 
interpretation of multiphase gas and kinematics in the 
simulation throughout this project, and Jessica K. Werk 
and Nicolas Lehner for helpful discussions on gas column 
densities. MEP acknowledges funding from NSF Awards 
AST-1008134 and AST-1312888, and the Clare Boothe 
Luce Program. MRJ acknowledges funding from AST- 
1312888. JEGP was supported by HST-HF-51295.01A, 
provided by NASA through a Hubble Fellowship grant 
from STScI, which is operated by AURA under NASA 
contract NAS5-26555. 
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